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Abstract. We present theoretical period - magnitude re- 
lationships for Cepheids in different filters for fundamental 
and first overtone pulsators, completing the work by Al- 
ibert et al. (1999). The results are provided for different 
metallicities characteristic of the Magellanic Clouds and 
the Milky Way. In contrast to the fundamental mode, we 
find a small metallicity effect on the period - luminosity 
relationship for the first overtone, due to the sensitivity 
of the period ratio P\/Po with metallicity. Comparison 
is made with observations from OGLE and EROS in the 
Small and Large Magellanic Clouds. We emphasize the 
constraint on theoretical predictions provided by the com- 
bination of both fundamental and first overtone observed 
sequences. We obtain excellent agreement between mod- 
els and data in a logP - Wj (Wesenheit index) diagram 
for a distance modulus for the LMC fio — 18.60 - 18.70. 
We analyse the uncertainties of the fundamental period - 
magnitude relationships and the consequences on distance 
determination. We show that an arbitrary shift of the in- 
stability strip by 350 K in T c s yields up to 0.45 mag effect 
on My at a given period, whereas the effect is less than 0.1 
mag in the if-band. Using recent near-IR observations in 
the Large Magellanic Cloud and our P - Mk relationship, 
we derive a distance modulus for the LMC in agreement 
with the value based on W\ data. 
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1. Introduction 

Over the last five years, a wealth of data for Cepheids 
was collected by several observational projects (HST key 
project, HIPPARCOS, microlensing experiments EROS, 
MACHO, OGLE, etc...). This fantastic amount of data al- 
lows unprecedented opportunities to understand the prop- 
erties of these variable stars. One of the relevance of 
Cepheids, and the main goal of the HST key project (cf. 
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Freedman 2000, and references therein), is the determina- 
tion of local extragalactic distance scales and consequently 
of the Hubble constant Hq. The period - luminosity rela- 
tionship (PL) of Cepheids is the conerstone of such pur- 
pose. An important contribution to the collection of data 
is provided by microlensing search projects toward the 
Magellanic Clouds (EROS: Renault et al. 1996; OGLE: 
Udalski, Kubiak and Szymanski 1997; MACHO: Alcock 
et al. 1997; and references therein). The recent release 
to the public domain of their Cepheid Catalog (EROS: 
Afonso et al. 2000; OGLE: Udalski et al. 1999a,b) provides 
large, homogeneous samples of high quality data which in 
principle allow a deep analysis of the period-luminosity 
relationships. In particular, a key question related to the 
PL relationship is its dependence on metallicity. This is of 
fundamental importance for extragalactic distance scales, 
usually based on an universal PL relationship. A contro- 
versy still exists regarding metallicity effects, both on the 
observational and theoretical viewpoints (see Alibert et 
al. 1999 and references therein). Another notorious de- 
bate concerns the distance to the LMC, which is an im- 
portant source of uncertainty in the determination of Hq 
(see Freedman 2000). Although PL relationships used for 
distance estimate are based on fundamental (F) pulsators, 
the afore-mentioned catalogs provide also statistically sig- 
nificant samples of first overtone (1H) Cepheids. Because 
of the distinction of F and 1H pulsator sequences in pe- 
riod - magnitude diagrams, 1H provides additional con- 
straints on theoretical models and on distance determina- 
tions based on PL relationships. 

In a recent paper, Alibert et al. (1999) performed self- 
consistent stellar evolution and linear stability analysis 
calculations for Cepheids. Their analysis was essentially 
devoted to fundamental mode pulsators and their period- 
magnitude-color relationships as a function of metallicity. 
The analysis of their results for first overtone pulsators 
and the comparison with observations were hampered at 
that time by the lack of significant samples of data. The 
aim of the present paper is to complete the work of Alib- 
ert et al. (1999) for first overtone pulsators and test their 
results against the observational constraints provided by 



2 



I. Baraffe and Y. Alibert: Period - magnitude relationships 



both F and 1H pulsators. In section 2, we derive PL re- 
lationships for 1H and in §3 we compare the results to 
observations. In §4, we briefly discuss the uncertainties of 
our period - magnitude relationships and consequences on 
distance determinations. 

2. Period - magnitude relationships for first 
overtone pulsators 

The description of the stellar evolution and pulsation cal- 
culations are given in details in Alibert et al. (1999). We 
recall the main ingredients: (i) evolutionary models for 
Cepheids are constructed with the Lyon evolutionary code 
from 3 to 12 Mq. We consider various initial compositions 
(Z,Y) = (0.02,0.28), (0.01,0.25) and (0.004,0.25) Q, rep- 
resentative of respectively the Galactic, the Large Mag- 
ellanic (LMC) and the Small Magellanic (SMC) clouds 
environments. Evolutionary calculations do not include 
core overshooting, (ii) A linear non-adiabatic stability 
analysis is performed on the complete evolutionary mod- 
els along the evolutionary tracks. This provides consis- 
tent mass- age-period-luminosity relations, (iii) Static at- 
mosphere models and their corresponding synthetic spec- 
tra are calculated for the same compositions used in (i), 
providing magnitudes and colors for a given (Z, M, L, 
T off ). 

2.1. Period - luminosity and Period - magnitude 
relationships in different filters. 

In order to derive statistical PL relationships, a mean po- 
sition is assigned to each mass, according to the time spent 
in different locations in the instability strip (IS) and a lin- 
ear least-square fit to these points is derived. As shown 
in Alibert et al. (1999), the mean position of a given stel- 
lar mass, accounting for its evolutionary time, is roughly 
located in the middle of the IS. 

The minimum masses m m i n undergoing a blue loop in 
the 1H instability strip are 3.25 Mq, 4 Mq and 5 Mq for 
respectively Z = 0.004, 0.01 and 0.02. In order to avoid 
biases due to the change of slope predicted near these 
minimum masses (see Alibert et al. 1999, their §3.4), we 
exclude in the present analysis m m \ n to derive the mean 
PL relationships. Because of the reduction of the blue loop 
extension toward lower masses, the slope of the PL rela- 
tionships for 1H can be affected if the fit is derived down 
to m m ; n , yielding steeper relations. Bauer et al. (1999) ob- 
served such a change of slope for fundamental pulsators 
in the SMC, but did not observe it for 1H, although pre- 
dicted by the models (see the discussion by Alibert et al. 
1999, §4.2.3). The relationships given in Tables and | are 
then derived from 3.5 Mq for Z = 0.004, 4.25 Mq for Z 
= 0.01 and 5.5 M© for Z = 0.02. Note that Alibert et al. 
(1999) derived preliminary PL relationships for 1H based 

1 Z is the metal mass fraction and Y the helium mass fraction 



Table 1. Coefficients of the log P (in days) - log L/Lq re- 
lationships (slope, zero-point) for first overtone pulsators 
as a function of metallicity. 



Z=0.02 (1.247, 2.659)) 
Z=0M (1.252, 2.646) 
Z=0.004 (1.257, 2.613) 



on the whole range of masses unstable in the IS (including 
^min) and obtained slightly different results than in Table 

0- 

The period - magnitude relationships in different filters 
are given in Table ^ for F and 1H modes. BVRI magni- 
tudes are based on the Johnson-Cousins system (Bessell 
1990) and JHK magnitudes are defined in the CIT sys- 
tem (Leggett 1992). Table || displays also the reddening- 
free Wesenheit index W\ = I — 1.55(V — J) (Madore and 
Freedman 1991), with the coefficient 1.55 (= Aj/E(V-I)) 
resulting from standard interstellar extinction curves {e.g., 
Schlegel et al. 1998). 

2.2. Effects of metallicity 

Figure displays the mean log Pi - log L/Lq relationships 
(Fig. |l|a) and the blue and red edges in a log Pi - T e g dia- 
gram (Fig. |l|b) as a function of Z. Excluding the minimum 
mass undergoing a blue loop in the IS, which decreases 
with Z, there is no noticeable effect of Z on the location 
and width of the IS in Fig. lb. The periods corresponding 
to m m i n in Fig. [l] are log Pi ~ 0.3, 0.1 and for respec- 
tively Z = 0.02, 0.01 and 0.004. We note however a small 
effect of Z on the PL relationship, with log L increasing 
by 0.04 for Z increasing from 0.004 to 0.02 at a given Pi. 

In the following, our analysis is restricted to period 
log' Pi > 0.3, since log Pi = 0.3 is the minimum period for 
Z=0.02 1H pulsators undergoing a blue loop (see Fig. 
An inspection of the relationships given in Table [l] shows 
that the slopes hardly depend on Z, but the zero points 
are affected by metallicity. Although small, the increase 
of L with Z translates in brighter magnitudes at a given 
Pi, with the largest effects in the VI J bands and the 
maximum effect in the / band. We note that My, M\ 
and Mj are brighter by 0.1 - 0.15 mag, Mh and Mk are 
brighter by 0.8 - 0.13 mag at a given Pi when Z increases 
from 0.004 to 0.02. Between Z=0.004 and Z = 0.01 the 
effect is ~ 0.1 mag in the VIJHK bands. Interestingly 
enough, the effect is the smallest in the B band, with a 
difference of Mb less than 0.1 mag between Z = 0.004 
and Z — 0.02. Opposite trends are found for fundamental 
mode PL relationships which show the largest Z effect in 
the B band, but small effects (less than 0.1 mag) at longer 
wavelengths (see Alibert et al. 1999, §5). 

This effect of Z on the PL relationship for 1H is related 
to the behavior of the period ratio Pi / Po as a function of 
Z. This ratio decreases as Z increases (cf. Baraffe et al. 
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Table 2. Coefficients of the log P (in days) - Magnitude relationships (slope, zero-point) for fundamental mode (upper 
row) and first overtone (lower row) as a function of metallicity. 



z 


M B 


My 


Mi 


Mj 


M„ 


M K 


Wi 


0.02 


(-2.658,-0.648) 
(-2.669,-1.561) 


(-2.905,-1.183) 
(-2.948,-1.990) 


(-3.102,-1.805) 
(-3.171,-2.526) 


(-3.256,-2.183) 
(-3.333,-2.844) 


(-3.346,-2.428) 
(-3.434,-3.042) 


(-3.367,-2.445) 
(-3.453,-3.057) 


(-3.407,-2.770) 
(-3.515,-3.355) 


0.01 


(-2.712,-0.655) 
(-2.699,-1.515) 


(-2.951,-1.153) 
(-2.972,-1.923) 


(-3.140,-1.769) 
(-3.192,-2.463) 


(-3.286,-2.157) 
(-3.356,-2.795) 


(-3.377,-2.411) 
(-3.461,-3.006) 


(-3.395,-2.428) 
(-3.481,-3.019) 


(-3.433,-2.725) 
(-3.534,-3.300) 


0.004 


(-2.707,-0.633) 
(-2.765,-1.430) 


(-2.939,-1.081) 
(-3.01,-1.807) 


(-3.124,-1.686) 
(-3.212,-2.341) 


(-3.262,-2.076) 
(-3.376,-2.678) 


(-3.351,-2.336) 
(-3.478,-2.898) 


(-3.369,-2.350) 
(-3.501,-2.906) 


(-3.411,-2.623) 
(-3.545,-3.170) 



1998): we find for Z = 0.004, P X /P Q ~ 0.73 and for Z = 
0.02, Px/Pq ~ 0.7. Since for a given L, Pq is essentially in- 
dependent of Z, Pi decreases because of decreasing Pi/Po 
ratio for increasing Z. Conversely, for a given Pi, L in- 
creases with Z. 

The behavior of P\/Pq with Z can be understood by 
analysing the regions in the star which contribute to the 
period of each mode. Such regions can be determined by 
the Epstein weight functions which are obtained by ex- 
pressing the eigen frequency of each mode by an integral 
over the whole star (Epstein 1950; Cox 1980). As already 
noted by Epstein (1950), the greatest contribution to the 
fundamental mode period is located between r/R ~ 0.7 
and r/R ~ 0.9. For the first overtone, two regions are im- 
portant: a first inner region located between r/R ~ 0.4 
and r/R ~ 0.7 and an outer zone with r/R > 0.85. These 
regions are characteristic of the mass range of interest (5- 
10 Mq). The inner region for 1H includes the opacity peak 
due to metals at T ~ 2 10 5 K and the outer region covers 
the H-Hc ionization zone. In the first region covering the 
metal opacity peak, an increase of the opacity k, related 
to an increase of the metallicity, yields an increase of T 
or a decrease of the density p at fixed r (or P) . Note that 
this region is dominantly radiative and the temperature 
gradient depends directly on Z since the radiative gradi- 
ent is proportional to k. Consequently, in this region, the 
sound speed C s oc (P/p) 1 / 2 oc T 1 / 2 is more sensitive to 
the metallicity than in other outer regions and increases 
with Z. The integral of the sound travel time r s = dr/C s 
through this zone is thus more sensitive to the metallic- 
ity, decreasing as Z increases, than in the other regions 
located at r/R > 0.7. The period ratio P\/Pq is related 
to the ratio of r s through the regions relevant for 1H to 
t s in the zones contributing to F. Because of the larger 
sensitivity to Z of t s for 1H, this ratio decreases as Z 
increases. We checked our arguments by calculating the 
integral of t s through the different regions of interest for 
different cases of M, L and T e g covering the instability 
strip and for different values of Z . 

The behavior of L as a function of Pi with Z explains 
the trends previously mentioned of the fluxes in different 
filters. In particular, it explains the low sensitivity of the 
log Pi - M-q relationship to Z because of the compensating 
effects of (i) higher L with increasing Z at a given Pi and 



(ii) lower P-flux for a given L and T e g as Z increases, 
because of more metallic line absorption in the B bandpass 
(see §5 of Alibert et al. 1999). 

Finally, because of the maximum effect in the I- 
bandpass, we note a non negligible effect of Z on the index 
Wi, which differs by 0.14 mag at long Pi up to 0.18 mag at 
shorter Pi between Z=0.004 and 0.02. When Z increases 
from 0.004 to 0.01, W\ gets brighter by - 0.13 mag. 

3. Period - magnitude diagrams: comparison with 
observations 

Comparison with observations are made with the EROS2 
and OGLE2 data for the MC. EROS 2 (Afonso et al. 2000) 
reports 239 1H Cepheids in the SMC and 113 in the LMC. 
OGLE2 provides ~ 800 1H pulsators in the SMC and ~ 
500 in the LMC. 

3.1. SMC Cepheids: VI -Bands 

Figure ^| displays the comparison between models with 
metallicity Z = 0.004 and the OGLE2 data. EROS 2 data 
are not shown for the sake of clarity, but are compara- 
ble to the OGLE2 data. We adopt the dereddened data 
provided by the OGLE2 catalog, with different redden- 
ing corrections for each fields (see Udalski et al. 1999b 
for details) and the same distance modulus p = 18.94 as 
used in Alibert et al. (1999). This value, suggested by 
Laney and Stobie (1994), yields general agreement be- 
tween models and various observations for fundamental 
pulsators (cf. Alibert et al. 1999). Unstable models during 
the first crossing phase (large open circles) and the core 
He burning blue loop (filled circles) are indicated (see Al- 
ibert et al. 1999 for details). As shown in Fig. ||, a gen- 
eral good agreement is found for 1H pulsators. The fun- 
damental mode sequence is indicated by the mean PL re- 
lationship, and compared to the relationship of Udalski et 
al. (1999c), recently revised by Udalski (2000). Note that 
these authors truncate their SMC sample at logP < 0.4, 
to avoid biases at short periods. The agreement between 
predicted and observed F relationships is thus excellent 
for logP > 0.4. As preliminary noted by Alibert et al. 
(1999) and confirmed by the data, observations for 1H 
do not extend above logP ~ 0.6, which corresponds to 
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Fig. 1. (a) Log P (in days) - L diagram for 1H pulsators 
and different metallicities. The symbols indicate the mean 
position of each mass in the IS as a function of Z . The 
mean PL relationships are plotted for Z =0.02 (solid line), 
Z=0.01 (dashed line) and Z=0.004 (long dashed-line). (b) 
Location of the blue and red edges of the 1H IS in a log P 
- T c ff diagram. The symbols correspond to the same Z as 
in (a). The minimum masses undergoing a blue loop in 
the IS are 3.25 M , 4 M Q and 5 M Q for respectively Z = 
0.004, 0.01 and 0.02. 



m ~ 5M Q with standard (no overshooting) evolutionary 
models. Note that models including overshooting would 
yield a lower mass for such period. Our linear stability 
analysis however finds unstable 1H modes (as well as F 
modes) for masses m > 5M©. Only non- linear calcu- 
lations can determine the dominant mode of pulsation. 
In the same vein, models predict that 1H pulsators with 
logP <J are in the first crossing phase. The significant 
number of such low period pulsators provides another 
strong constraint on non-linear calculations which could 
test the following scenarios: (1) models on the first cross- 
ing, characterised by a lower luminosity than during the 
He core burning blue loop for a given mass, should favor 
1H as the dominant mode, and (2) models with m > 5M Q 
on the blue loop should oscillate predominantly in the fun- 
damental mode. 

We note however that in case (1), one may expect a 
change in the number density of objects with log P < 0, 
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Fig. 2. Period (in days) - Magnitude diagrams for 1H pul- 
sators in the VI bands for models with Z=0.004 and SMC 
observed Cepheids. The filled circles correspond to unsta- 
ble 1H modes during core He burning phase from m m i n 
= 3.25 Mq to 6 Mq. The large open circles correspond 
to first crossing unstable 1H modes for 3, 3.25, 3.5 and 4 
Mq. Observations (dots) are from OGLE2 (Udalski et al. 
1999b). The distance modulus is 18.94. The 2 solid curves 
are the theoretical mean PL relationships for respectively 
F (lower curve) and 1H (upper curve) pulsators. The long- 
dashed curves are the observed PL relationship derived by 
Udalski et al. (1999c) and Udalski (2000) for F and 1H 
pulsators. 

given the much faster evolutionary timescale on the first 
crossing compared to the blue loop phase. Such a change 
is not displayed by current observations. A detailed statis- 
tical analysis taking into account evolutionary timescales, 
mass functions and eventual observational biases is re- 
quired to investigate this point. 

3.2. LMC Cepheids: VI -Bands 

Comparison between models with Z=0.01 and observa- 
tions in the LMC is shown in Fig. |3| Data are from EROS2 
(Afonso et al. 2000) and OGLE2 (Udalski et al. 1999a). 
We adopt the same reddening corrections as in Udalski et 
al. (1999a) for the OGLE2 sample, which varies from field 
to field. For the EROS2 data, we adopt the same redden- 
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ing correction E(B — V)=0.10 as used in Alibert et al. 
(1999). We note that adopting a constant E(B - F)=0.10 
for the OGLE2 data has no noticeable effect when compar- 
ing models and data in Fig. |[ We adopt the same distance 
modulus 18.50 as in Alibert et al. (1999). The theoretical 
PL relations for F and 1H modes are compared to the ob- 
served relations derived by Udalski et al. (1999c) on the 
LMC sample limited to logP > 0.4, and recently revised 
by Udalsky (2000). The same agreement as for the SMC 
data is found and the models reproduce the width and 
location of the instability strip satisfactorily. We however 
note discrepancies between models and observations for 
the mean P - My relationship for both F and 1H modes. 
This discrepancy appears also, but to a lesser extent, for 
the SMC (cf. Fig. |) and is better illustrated in P - Wj 
diagrams (see next section). 

As illustrated in Fig. 3, no 1H pulsators are observed 
above logP ~ 0.8, corresponding torn > 7M . Note that 
the lack of 1H pulsators at long periods could also be in- 
terpreted in terms of very small amplitudes below the level 
of detection of OGLE or EROS. However, there is no ob- 
vious decrease of the amplitudes of 1H pulsators in the 
LMC, as well as in the SMC (see Afonso et al. 1999), as 
the period increases which could support this interpreta- 
tion. Below logP ~ 0.1, observations can be explained 
by first crossing models. We therefore find the same type 
of constraints as derived in the previous section from the 
SMC data on non-linear calculations for the determina- 
tion of the dominant pulsating mode. Such comparisons 
suggest that the afore-mentioned properties in scenarios 
(1) and (2) (see §3.1) are intrinsic to 1H pulsators for Z 
varying from 0.004 to 0.01. More observations at higher 
metallicities are currently required to determine if these 
properties apply also to Galactic Cephcids. 

3.3. P - W\ diagram 

Figures |J and [| display observations and models for both 
F and 1H in a log P - Wj diagram, for respectively the 
SMC and the LMC. The Wesenheit index W\ is an useful 
quantity since it is rather insensitive to the reddening and 
can in principle remove part of the scatter due to differ- 
ential reddening. For the LMC, the Z — 0.01 models are 
also compared to observations from Gieren et al. (1998) 
for F Cepheids. Note that Alibert et al. (1999) found good 
agreement in different optical and near-infrared P - mag- 
nitude diagrams with the latter sample of data. As ex- 
pected, the scatter of data in such diagram is smaller than 
in the V and / bands. A first inspection of Figs. |I] and || 
show a general agreement between models and observa- 
tions for both SMC and LMC: the models corresponding 
to different masses cover the observed location of both F 
and 1H Cepheids. We note that the models are in better 
agreement with the observed width of the IS in the LMC, 
whereas the data in the SMC show a larger scatter. This is 
consistent with the recent analysis by Groenewegen (2000) 
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Fig. 3. Period (in days) - Magnitude diagrams for 1H pul- 
sators in the VI bands for models with Z=0.01 and LMC 
observed Cepheids. The filled circles correspond to unsta- 
ble 1H modes during core He burning phase from m m i n — 
4 Mq to 7 Mq . The open circles correspond to first cross- 
ing unstable 1H modes for 3.875, 4, 4.25, 4.5, 5 M Q . Ob- 
servations (dots) are from OGLE2 (Udalski et al. 1999a) 
and EROS 2 (Afonso et al. 2000). The distance modulus 
is 18.50. The 2 solid curves are the theoretical mean PL 
relationships for respectively F (lower curve) and 1H (up- 
per curve) . The 2 long-dashed curves are the observed PL 
relationships derived by Udalski et al. (1999c) and Udalski 
(2000). 



who interprets the larger dispersion in the PL relation in 
W\ for the SMC compared to the LMC in terms of a larger 
intrinsic depth of the former Cloud. 

For both the SMC and the LMC, the predicted log P 
- W\ mean relationships are shifted compared to the ob- 
served relationships. Since the shift between observed and 
predicted relationships is almost constant on the whole 
range of P and is the same between F and 1H relations, 
we find that an increase of the distance moduli fio by the 
same amount ~ 0.15 - 0.2 mag for both clouds removes 
easily this discrepancy. This is illustrated in the insets 
of Figs. ^ and ||. Note that such an increase of does 
not alter the agreement found between observations and 
models in the /-band (see Fig. 2 and Fig. 3). It however 
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increases the discrepancy in the U-band mentioned in the 
previous section. It is difficult to analyse the reason of 
such inconsistency. The choice of the reddening and of the 
extinction curve coefficients, affecting mostly the U-band, 
could be the reason of this problem. This emphasizes the 
large uncertainties intrinsic to observations in this filter 
and suggests to use rather W\ to derive a distance mod- 
ulus. Note however that, independently of reddening and 
distance modulus, the models do predict the correct lo- 
cation of the observed 1H sequence relative to the F se- 
quence, for both SMC and LMC. 

We also note that the models yield a difference of 
distance moduli between SMC and LMC of ~ 0.45, in 
good agreement with previous determinations based on 
Cepheids (Laney and Stobie 1994; Udalski et al. 1999c), 
RR Lyr or Red Clump stars (cf. Udalski et al. 1999c; Udal- 
ski 2000). Finally, Caputo et al. (2000) also derived theo- 
retical log P -W\ relationships for the fundamental mode 
which are shallower than the present relationships. Al- 
ibert et al. (1999) already mentioned and discussed this 
difference for fundamental mode period - magnitude rela- 
tionships. As illustrated in Figs. [| and [|, the Caputo et 
al. (2000) log P - W\ relationships are also shallower com- 
pared to the observed data in both SMC and LMC. For the 
distance moduli adopted in the insets of Figs. [| and ||, the 
F sequence of Caputo et al. (2000) reproduce the F data 
at long period, but merges into the observed 1H pulsators 
at shorter periods. A variation of the distance modulus 
cannot clearly solve such discrepancy. This highlights the 
constraint on models and distance moduli provided by the 
combination of F and 1H observed sequences. 

4. Uncertainties on distance determination 

Alibert et al. (1999) analysed the effects of uncertainties 
inherent to stellar evolution models: convection treatment, 
overshooting, mass loss, initial helium abundance, etc... 
The main conclusion is that such uncertainties barely af- 
fect the PL relationships and the conclusion on mctallicity 
effects. We however recall that such uncertainties can af- 
fect the comparison between models and observations if 
the mass is used as another constraint. 

The main uncertainty in our calculations is due to the 
neglect of convection - pulsation coupling in our linear sta- 
bility analysis. Convection is frozen in, which means that 
the perturbation SF conv of the convective flux is neglected 
in the linearized energy equation. An arbitrary criterion 
is used to define a red edge for the instability strip (see 
Alibert et al. 1999 for details), since such approximation 
cannot yield naturally a red edge. 

A time dependent non local theory of convection is 
required to take into account the effect of convection on 
pulsation. Such a theory is lacking and the current recipes 
include several free parameters (cf. Yecko et al. 1998 and 
references therein). Despite the lack of a robust theory, it 
has now become clear that convective energy transport is 



essential to describe the observed properties of pulsating 
stars such as light curves or Fourier coefficients (Yecko et 
al. 1998; Feuchtinger 1999; Feuchtinger et al. 2000). 

In a systematic analysis based on linear stability anal- 
ysis of Cepheid models, Yecko et al. (1998) have shown the 
high sensitivity of the position of the IS to these free pa- 
rameters, which need to be calibrated according to astro- 
nomical observations. But as emphasized by Feuchtinger 
et al. (2000) even a correct description of the instability 
strip requires hydrodynamical calculations, since the red 
edge is determined by nonlinear effects. Even for the blue 
edge, Feuchtinger et al. (2000) report a shift toward higher 
T c ff by about 350 K for the fundamental mode when con- 
vection is included in the pulsation calculation. 

In the present paper, we derive a rough estimate of 
the sensitivity of the period - magnitude relationships de- 
rived under the afore-mentioned approximations, and the 
effects on distance determination. A more detailed analy- 
sis and the quantification of uncertainties resulting from 
convection is addressed in a forthcoming paper (Alibert & 
Baraffe 2001). Since our results yield a general good agree- 
ment with the width of observed instability strips (see Al- 
ibert et al. 1999), a first estimate of such uncertainties 
can easily be derived by shifting arbitrarily both the blue 
and red edges of the instability strip toward cooler or hot- 
ter T c ff, keeping the width unmodified. In the following, 
we focus on fundamental PL relationships, since they are 
mostly used for distance determinations. Inspired by the 
afore- mentioned results by Feuchtinger et al. (2000), we 
adopt a shift in T e g of 350 K. Based on the models appro- 
priate for the LMC (Z = 0.01), a shift by 350 K toward 
hotter (cooler) T e g yields brighter (fainter) magnitudes at 
a given P. The effect decreases toward ncar-IR filters: the 
resulting variation of the magnitude in the i?-band reaches 
up to 0.6 mag, it does not exceed 0.3 mag in the /-band 
and remains below 0.2 mag at longer wavelengths. This 
is expected, given the decreasing sensitivity to T e g of the 
flux toward longer wavelengths. In the U-band, My varies 
from 0.3 mag at short P (logP = 0.5) up to 0.45 mag at 
longer P (logP = 1.7). In the K band, the effect is less 
than 0.1 mag on the whole range of periods. This simple 
test highlights the high sensitivity of distance determina- 
tions based on theoretical P - My relationships, since a 
variation of 350 K on the location of the blue or red edge 
can easily result from uncertainties due to convection (see 
Yecko et al. 1999; Feuchtinger et al. 2000; Alibert and 
Baraffe 2001). 

Finally, since our simple test clearly shows that near- 
IR period - magnitude relationships are more reliable, it 
is worth to derive a distance modulus for LMC based on 
near-IR data. Recent observations in the near-IR (Madore, 
2000 priv. comm.; Persson et al. 2000) complete data al- 
ready available from Laney and Stobie (1994) and Gieren 
et al. (1998). The P - M K relationship (Table 2) yields a 
distance modulus for LMC of 18.60- 18.70, as illustrated 
in Fig. Interestingly enough, this distance modulus is in 
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Fig. 4. log P (in days) - W\ diagram for F and 1H pulsators for models with if=0.004 and SMC observed Cepheids. 
Symbols and curves are the same as in Fig. 0. The inset displays only the mean relationships with the same linestyle 
as in the main figure for a distance modulus fio — 19.14. The dash-dotted curve corresponds to the F relationship of 
Caputo et al. (2000). Note that for 1H, the theoretical and observed mean relationships are indistinguishable in the 
inset. 



good agreement with the value determined in §3.3 based 
on W\ and the combination of F and 1H observed se- 
quences. 

5. Conclusion 

We have derived period - magnitude relationships for 1H 
pulsators in different filters. In contrast to fundamental 
mode, we find a small effect of metallicity in PL relation- 
ships for first overtone pulsators. This effect is due to the 
dependence on metallicity of the ratio P\/Pq. Our mod- 
els reproduce the location of both F and 1H observed se- 
quences of OGLE2 in the SMC and the LMC. Using the 
reddening free index W\, models and observations are in 
good agreement for a LMC distance modulus fio = 18.65 
- 18.70. We note however an inconsistency since this value 
yields a significant discrepancy between predictions and 
data in the V-band. Such inconsistency may illustrate 
problems in the choice of the extinction coefficients and 
the reddening, which are important for the comparison 



between models and data in the U-band but not crucial 
for Wi. 

We show that an arbitrary shift of 350 K in T c g of the 
location of the instability strip yields up to 0.45 mag effect 
on My at a given P. Toward near-infrared wavelengths, 
the effect is smaller and less than 0.1 mag in the if -band. 
Using recent Cepheid data in the K-band for fundamen- 
tal pulsators in the LMC, the models predict a distance 
modulus for LMC /x = 18.60 - 18.70, in good agreement 
with the predictions based on Wi and combined 1H and 
F data. 

Although Cepheid observations in the U-band are 
widely used for distance determination, we stress the high 
uncertainties inherent to observations in this filter (red- 
dening correction, intrinsic dispersion in the instability 
strip, etc..) and the high sensitivity of P - My relation- 
ships to theoretical uncertainties. We also emphasize that 
even if the coupling between convection and pulsation pro- 
vides a large source of uncertainty in the present work 
and in all current theoretical calculations, predictions in 
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Fig. 5. log P (in days) - W\ diagram for F and 1H pulsators for models with Z=0.01 and LMC observed Cepheids. 
Symbols and curves are the same as in Fig. ||. The triangles are observations from Gieren et al. (1998). The inset 
displays only the mean relationships with the same linestyle as in the main figure for a distance modulus = 18.70. 
The dash-dotted curve corresponds to the F relationship of Caputo et al. (2000). 



the X-band are expected to be much less affected by such 
uncertainty. 
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